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Abstract. Infrared spectra of two type II supernovae 6 months after explosion are presented. The spectra exhibit 
a strong similarity to the observations of SN 1987A and other type II SNe at comparable epochs. The continuum 
can be fitted with a cool black body and the hydrogen lines have emissivities that are approximately those of a 
Case B recombination spectrum. The data extend far enough into the thermal region to detect emission by the 
first overtone of carbon monoxide. The molecular emission is modeled and compared with that in the spectra of 
SN 1987A. It is found that the flux in the CO first overtone is comparable to that found in SN 1987A. We argue 
that Carbon Monoxide forms in the ejecta of all type II SNe during the first year after explosion. 
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The presence of molecular emission in the spectra of super- 
novae while not new remains rarely detected. First detec- 
tion of carbon monoxide in a supernova was in the spec- 
tra of S N198 7A (Catchpole & Glass |1987| McGregor & 
Hyland 1987, Ames Research Center 1987, Oliva et al. 



1987, Spyromilio et al. 1988). In SN 1987A observations 



of the first overtone of CO at 2.3-/im were complemented 
by observations at 4.6-/im of the fundamental bands. In 
addition to carbon monoxide, the spectra of SN 1987A re- 



vealed bands of SiO (Meikle et al. |1993| ; Roche et al. [1991]) 
and CS (Meikle et al. |1993| ). The presence of H^ i s also 
claimed in the spectra of SN 1987A (Miller et al. |1992[) 



although alternative atomic species can also explain the 



identified features. Spyromilio & Leibundgut (1996) re- 



ported the detection of CO first overtone emission in the 
spectra of the type II supernova 1995ad while Gerardy et 
al. (2000) and Fassia et al. ( |2000| ) reported the presence 
of CO in SN 1998S. 

The formation of molecules in the ejecta of supernovae, 
even CO - the most stable of diatomic molecules with a 
dissociation energy of 11.09eV (Douglas & M0ller 1955Q , 
is not trivial. The high UV field due to recombinations 
and the energetic electrons from the radioactive decays of 
56 Ni and its daughter 56 Co create an inhospitable environ- 
ment with multiple dissociation paths (see Lepp, Dalgarno 



& McCray 1990). The presence therefore of molecules in 
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the ejecta imply a particular distribution of the material 
within the ejecta and place constraints on the degree of 
mixing possible. These effects have only been analyzed in 
detail in the well observed SN 1987A. The degree to which 
that object is representative of the entire class of type II 
SNe, which will almost by definition not be observed with 
the same accuracy, is of some importance. 

Molecules also affect the degree of cooling of the ejecta. 
The partition function of a diatomic molecule such as CO 
is enormous compared to that of even a heavy ion such 
as Fe + . Moreover, the lower rotational levels of the fun- 
damental, first and second overtones of CO require low 
energies for their excitations allowing them to continue 
cooling the ejecta even as the temperature drops. 

Gerardy et al. ( |200C ) and Fassia et al. ( 200C| ) also re- 
port the detection of emission from dust in SN 1998S and 
Fassia et al. argue that the cooling by CO may lead to 
conditions that favour the formation of dust. 

Here we report on observations of two more type II 
SNe in the near infrared which exhibit emission by carbon 
monoxide. In sect. 2 we report on the observations. The 
data are discussed and compared with those of SN 1987A 
and theoretical models in sect. 3. 



2. Observations 

SN1998dl in NGC1084 (reces sion velocity 1411 kms" 1 ) 
was discovered by King et al. ( 1998 ), on 1998 August 20, 
from observations made at the 0.8-m Katzman Automatic 
Imaging Telescope (KAIT). It was subsequently found to 
have been present in i mages taken on 1998, August 2 and 
August 9. Filippenko ( 1998; ) detected hydrogen Balmer 
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lines in early optical spectra classifying the object as a 
type II supernova. 

SN1998dl was observed at the ESO 3.58-m New 
Technology Telescope on La Silla, on 1999 January 1, at 
an age of approximately 150 days. The near-infrared cam- 

|199SD 



was 



era and spectrograph SofI (Moorwood et al 
used for these observations. The red low-resolution grism 
was used for the H and K band spectra. While the grism 
does cover the region between 1.8 and 1.92 /^m the atmo- 
spheric transparency is very low and no useful data are 
available in this region. A 1 arcsecond slit was used which 
gives a resolution of A/AA of 650 and a wavelength cov- 
erage from 1.5 /im to 2.5-/im. The spectroscopic standard 
BS 2290 (Allen & Cragg 1983| ) was used to remove atmo- 
spheric features and to correct for the spectral response of 
the instrument. The conditions were not photometric and 
therefore the absolute fluxing of the spectrum is arbitrary. 
The slit projects to 3 pixels on the detector oversampling 
the data. Narrow features, of order 3 pixels (e.g. the par- 
ticularly prominent feature at 2 /im) , in the spectra are 
not real but result from incomplete corrections for atmo- 
spheric features. 

Supernova 1999em in NGC 1637 (recession velocity 
717kms _1 ) was discovered by Li ( 1999 ) from observ ations 
made, on 1999 Oc tober 29, at the KAIT. Jha et al. ( |l999D 
and Deng et al. ( |1999| ) classified the object as a type II 
supernova based on spectroscopic observations. 

SN 1999em was observed, on 2000 April 24, at an age of 
approximately 170 days with the same instrumental con- 
figuration as 1998dl. The spectroscopic st andar d used in 
this case was HIP 27855 (Perryman et al. 1997). The ac- 



curate removal of atmospheric features requires that the 
standard and the target are observed using the same re- 
solving power. The usage of a narrow slit in both observa- 
tions of the standard and the target implies that a sepa- 
rate photometric calibration of the data needs to be made 
to recover absolute fluxes. The imaging mode of SofI was 
used for this p urpos e and the NICMOS standard S121-E 
(Persson et al. 1998 ) was used for the calibration. The de- 
rived K s magnitude of the supernova was 15.1. The mean 
flux of the observed spectrum was corrected to the this 
magnitude. 



3. Discussion 

The spectra of SN 1998dl and SN 1999em are displayed in 
Figs. G] and 0, respectively. The data are shown at their 
observed wavelengths without correction for the recession 
velocity of the parent galaxies. The spectra are character- 
ized by a strong continuum upon which emission lines are 
superimposed. The strong Brackett7 line (rest wavelength 
2.1655/im) is evident. In the SN 1999em data the shorter 
wavelength transitions of the Brackett series can also be 
seen. Long-wards of the Br7 line, at 2.2 1 /im, a broad fea- 
ture identified by Meikle et al. ( 1989 ) in the spectra of 
SN1987A as emission by the 2.207 /mi Nal 4s 2 S-4p 2 P° 



multiplet is evidentn. Starting at 2.3 /mi and extending to 
the end of our spectral converage a broad feature is evi- 
dent in the data. This broad feature we attribute to carbon 
monoxide first overtone transitions (Av = 2). The much 
weaker second overtone (Av — 3) is expected to start at 
1.5 /an. Its contribution is expected to be only 1.4% of 
the strength of the first overtone (Bouanich & Brodbeck 



1974) and is therefore not detectable in our spectra. In 



much higher signal to noise data as would possibly be ob- 
tained at the VLT with ISAAC it might be possible to 
detect such emission. The fundamental (Av = 1) band 
at 4.6 /im is very much stronger than the first overtone 
observed here but also in an extremely unfavourable at- 
mospheric window. This fundamental band, as mentioned 
in the introduction, almost dominates the cooling of the 
ejecta for some time and observation would be of great in- 
terest. Unfortunately they are probably limited to a com- 
bination of the very biggest telescopes and the very nearest 
supernovae. 

A strong emission feature is observed at 2 /im. As noted 
above the narrow spike superimposed on the broad emis- 
sion is instrumental. In SN 1987A, a broad feature at the 
same wavelength was identified by Meikle et al. as emission 
by the Cal 4p 3 P°-3d 3 D multiplet and the [Fel] a 5 F 5 -a 3 F 4 
multiplet. Assuming a contribution to the 2/tm by [Fel] 
in the spectrum of SN 1998dl, the strong line at 1.51 /tm 
could be identified as emission by the [Fel] a 5 D— 2-a 5 F4 
multiplet. Alternatively Mgl also has a strong transition 
at 1.5031 /tm. The lower recession velocity of the parent 
galaxy of SN 1999em prevents us from detecting the same 
line in that spectrum. 

At 2.06 /im both in SN 1998dl and SN 1999em a small 
absorption is present. This we identify as the P-Cygni 
trough of the Hel 2s 1 S-2p 1 P° 2.058 /tm transition. We 
note that this region is one of poor atmospheric transmis- 
sion with a strong telluric CO2 feature at 2.06 /tm. While 
we believe that we have correctly compensated for this 
it is possible that the feature we identify as due to Hel is 
due to poor cancellation by of aforementioned atmospheric 
feature. Hel is expected to be very weak but its presence 
would be significant as highly excited He is a tracer for 
the presence of energetic electrons from the radioactive 



decay (Graham |1988| Fassia & Meikle |199J). The singlet 
transition observed here would need to develop sufficient 
optical depth to exhibit a P-Cygni trough. While transi- 
tions from 2s 1 S to the ground state of Hel are forbidden, 
depopulation of this energy level can occur via two photon 
decay with a low transition probability (A^ls^ 1 ). Given 
the age of the supernovae at the time of observation even 
this process would have de-populated the ground state of 
the 2.058 /tm line. Recombinations are the obvious source 
of electrons which, as the transitions from the 1 P° states 
to the ground are expected to be saturated, will naturally 
cascade down to the 2s 1 S level. Given the high ionization 
potential of He (24.587 eV) the presence of the 2.058 /tm 



1 In Meikle et al. ( L989J) the lower level of this multiplet in 
table 3 is incorrectly labeled as 4s 2 D 
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Fig. 1. Spectrum of 1998dl. 
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Fig. 2. Spectrum of 1999em. 



transition suggests that the excitation by energetic elec- 
trons argued by Graham and Fassia & Meikle is occuring 
also in these two SNe. 

We have used simple line models to make fits to the 
data. These fits are shown in Figs. [| and [|. The ingredients 
of the models are described below. 

The continuum in the models is a black-body with 
a temperature of 3300 K for SN 1998dl and 2000 K 
for SN 1999em. The flux observed in the spectrum of 
SN 1999em at this temperature corresponds to a surface 
of last scattering expanding at 500 kins -1 for 170 days at 



a distance of 7.8 Mpc (Sohn & Davidge |1998|) . 

For the carbon monoxide a simple LTE (Boltzmann 
distribution) model of the populations of the rotational- 
vibrational energy levels in the ground state electronic 
level x £ of 12 C 16 (CO) was used. T his mo del has been 
described before in Spyromilio et al. ( 1988 ). Briefly, the 



energy levels for a given rotational level J and a given 
vibrational level v can be determined using: 



E = lu(v + i) - lux(v + i) 2 + uy(v + ±) 3 
+B v J(J+l)-D v J 2 (J+lf 



(1) 



with Lo 1 ujx,ujy,B v and D v from Young (1968). The A 
values are given by 

64ttV 



A , — 



3hgj 



-{tmy 



(2) 



where v is the transition frequency, gj = 2 J + 1 and 
TM is the transition matrix from vibrational level v to v' . 



TM for AJ of ±1 is given by Chackerian & Tipping fll983J) 
The lowest 15 vibrational levels were included and within 
each 110 rotational level populations were calculated. 

As has been sho wn in the spectra of SN 1987A 
(Spyromilio et al. 1988J) the CO band heads are sensitive to 
the excitation temperature. Higher temperature in the CO 
model increase the contribution of higher vibrational pairs 
and increase the emission long- wards of the v = 2 to v' = 
band-head. The quality of the data is not good enough 
to reliably distinguish between small changes in the tem- 
perature. The SN 1998dl data indicate a slightly higher 
temperature than the SN 1999em. The fits use 2500 K for 
SN1998dl and 2000 K for SN 1999em. The uncertainty in 
these values is of order a few hundred degrees K. The dif- 
ference in the temperatures leads to a significantly differ- 
ent appearance of the first overtone emission. The 2000 K 
temperature for the SN 1999em spectrum is based on the 
assumption that the drop in flux long-wards of 2.35 //m 
is real. The SN 1998dl spectra appear to show a much 
smoother decline towards longer wavelengths. 

The velocity at which the CO is expanding can be used 
to co nstrain the mass of the progenitor (Gerardy et al. 
|2000 ). In SN 1998S Gerardy et al. derive a high expansion 
velocity (v ~ 3000 kins -1 ) based on the smooth rise of the 
2-0 R-branch band h ead an d a progenitor mass in excess of 
25 Mq. Fassia et al. (200C) do not detect such a steep rise 
and argue for a lower velocity. In addition their modeling 
of the ejecta favours a lower mass progenitor. We argue 
below in favour of the low range of velocities although we, 
like Gerardy et al., have a smooth rise of the 2-0 R-branch 
band head. 

In the spectra of SN 1987A a rigorous fitting of the ex- 
pansion velocity could be made (Spyromilio et al. 1988) 



based on the appearance of the R-branch emission from 
each pair of e nergy states. Both Gerardy et al. ( [2000 ) and 
Fassia et al. (2000) use the lack of clear separation of the 
R-branch band heads to place lower limits on the expan- 
sion velocity in the observations of SN 1998S. This same 
effect can be used by us to place a weak lower limit of 
1000 kms -1 on the expansion velocity. 

The comparison of SN1998S with SN 1987A (see 
Gerardy et al. Fig. 6) shows that another line is possi- 
bly blended with the CO band head at 2.28 /jm. This fea- 
ture, which at the recession velocity of SN 1987A is ob- 
served at 2.26 fim, is discussed later but remains uniden- 
tified. The 2.26 /im feature is very strong in the spectra of 
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SN 1999em and weakly present in our data on SN 1998dl. 
The same feature was present in the spectra of SN 1995ad 
(Spyromilio & Leibundgut 1996). We cannot account for 
the lack of detection of this in the data of Fassia et al. 
(200C). The blending of this feature into the band head 



SN 1998cll SOFI I ITT 999-01-02 



implies that models that do not take this into account 
will underestimate the true sharpness of the band head 
rise and derive higher than expected expansion velocities. 
The best estimate of the velocity of the CO, we believe, 
comes from the comparison of the observed spectra with 
those of SN 1987A (see Figs. | and |). The good matching 
of the rise of the 2-0 band head suggests that the veloc- 
ity o bserve d in SN 1987A (1800-2000 kms" 1 ; Spyromilio 
et al. 1988 ) is representative of the velocity of the CO in 
SN 1998dl and SN 1999em. Clearly in the absence of the 
blended line we would agree with Gerardy et al. ( 2000 ) 
and also derive a much higher expansion velocity for the 
CO. Each CO line in our models was convolved with a 
Gaussian line profile of FWHM 2000 kms" 1 . 

The Brackett series line strengths in the model as- 
sume a case-B recombination spectrum at a temperature 
of 3000 K and electron density of 10 5 (Hummer & Storey 
1987). The Br7 and 9^4 transitions are well fit using 



these parameters. The BrS line lies at the very edge of 
the atmospheric window and the cancellation of the at- 
mospheric features there is poor. While the norm would 
be not to display the spectra in this region due to the 
large uncertainties in the data we include them here. The 
data show that the spectra are consistent with the model. 
We do not derive any further information from the region 
around 1.92 mum. Unfortunately the Bre transition lies in 
the gap between the two atmospheric windows and is not 
observable from the ground. The FWHM of the Gaussian 
line profile computed was 3000 kms -1 . As has been the 
case in the other SNe in which CO was detected, the ve- 
locity of the hydrogen lines is higher than that of the CO. 

The individual lines identified earlier as Nal, Cal and 
[Fel] have not been included in the displayed fits as the 
individual contributions of the species would be arbitrary. 

As already discussed with respect to the expansion ve- 
locity, the observations presented here can be compared 
wi th the observations of SN 1987A (data from Meikle et 
al. |1989D . Fig. | shows the data from SN 1987A taken at an 
age of 192 days superimposed on the spectra of SN 1998dl. 
The SN 1987A data have been scaled to match the flux in 
the CO feature and red-shifted to the recession velocity of 
NGC1084. The continuum underlying the SN 1987A data 
has been artificially adjusted to match the average con- 
tinuum level in the SN 1998dl data. The SN 1998dl con- 
tinuum is much bluer than that of SN 1987A at the epoch 
at which the comparison data are taken. The poor correla- 
tion at 2.25 (im is partially due to this effect and partially 
due to the difference in the strength of the 2.26 fim feature. 
Given the quality of the SN 1998dl data, this comparison, 
adds plausibility to the identification of CO in the spec- 
tra. Also as mentioned above it allows a more accurate 
determination of the velocity structure of CO than can be 
achieved from the models. A similar comparison is made 
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Fig. 3. Spectrum of 1998dl and model. Features discussed 
in the text are identified in the figure. 
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Fig. 4. Spectrum of 1999em and model. Features dis- 
cussed in the text are identified in the figure. 



in Fig.| for SN 1999em using SN 1987A data obtained at 
an age of 255 days. 

The flux within the CO emission band for SN 1999em 
is 5.5xl0~ n ergs/s/cm 2 which at th e distance of NGC 
1637 (7.8 Mpc; Sohn & Davidge |l99^ ) is 20% lower than 
that of SN 1987A at a comparable epoch. 

The mass that corresponds to this emission depends 
critically on the emission model. Our LTE model gives 
masses of order 10 -4 M Q . Liu & Dalgarno 1995 predict a 
much higher mass (^1O _3 M0) for the same emission at 
this age and argue that the conditions within the ejecta 
favour the formation of this amount of CO. In fact it is 
expected that significantly more CO formed during the 
youth of the supernova (« 100 days after explosion ) but 



has been destroyed. Gearhart, Wheeler & Swartz 1999 



favour the formation of much less CO than Liu & Dalgarno 
and more in line with our mass estimates. Both Gerardy 
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Fig. 5. Spectrum of 1998dl in the K band with spectrum 
of SN1987A superimposed. 

SN 1999em SOFI NTT 2000-04-24 & SN 1 987A at 255 days 
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Fig. 6. Spectrum of 1999em in the K band with spectrum 
of SN1987A superimposed. Note that the feature marked 
as CO + is much weaker in 87A than in 99em. 



et al. fl2000D and Fassia et al. (|2000|) employ non-LTE and 
optically thick models to estimate the mass of CO emitting 
and derive high masses. These models base their input on 
assumptions and parameters that cannot be derived from 
our limited observations and even in the case of SN 1998S 
are based on assumed filling factors. Moreover Gerardy et 
al. argue that their comparisons are more a proof of prin- 
ciple rather than true models. We restrict ourselves there- 
fore to our comparison with SN 1987A. Complications due 
to clumping and asphericity are beyond the scope of this 
work to address but it is worth noting that polarization 
has been detected in the spectra of SN 1999cm (Leonard, 
Filippenko & Chornock tL999t Wang et al. 119991). 



Sollcrman et al. ( 2001 ) detected no sign of dust in spec- 
troscopic observations of SN 1999em obtained at an age of 
450 days. While the formation of molecules may provide 



the cold sites for dust formation in the ejecta of type II su- 
pernovae and therefore may be a necessary precondition, 
it is clearly not sufficient. 

The feature at 2.26-/xm short-wards of the CO v = 2 — > 
band head remains unidentified. The formation of suffi- 
cient CO + , whose first overtone v = 2 — > band head co- 
incides with this feature, is excluded by both Gearhart et 
al. and Liu & Dalgarno. Note, however, that Petuchowski 

al. " 



et 



1989 



have argued that CO + could form in suffi- 
cient quantities. The feature was reported by and iden- 
tified as CO + by Spyromilio et al. (1988) and reported 
by Spyromilio & Leibundgut in the sp ectra of SN 1995ad. 
It is also obse rved by Gerardy et al. fl2000|) in SN 1998S. 
Fassia et al. ( 2000 ) did not detect the emission in their 
spectra of SN 1998S. While theoretical models fail to pre- 
dict the necessary amounts of CO + the identification of 
this feature remains insecure. 

4. Conclusions 

The spectra presented in this paper combined with the ob- 
servations presented by Spyromilio & Leibundgut (1996), 
Gerardy et al. ( 200C ) and Fassia et al. ( 200C ) suggest that 
the formation of CO in the ejecta of type II supernovae is 
ubiquitous at an age of 3 to 6 months. 

The mass of CO present in the ejecta is difficult to es- 
timate accurately. Our LTE models while simplistic rely 
on few assumptions with respect to the distribution and 
excitation of the CO. While for well studied supernovae 
such as SN 1987A the much more detailed models of Lepp, 



Dalgarno and McCray (1990) are well justified, it is ar- 
guable whether the input parameters for such models are 
well enough constrained in the case of the objects de- 
scribed here. We therefore argue that since the absolute 
fluxes argee for SN 1987A and SN 1999em so, to first ap- 
proximation, will the emitting masses. 

The expansion velocity of the CO is not accurately 
determined. By comparison with SN 1987A we argue that 
velocities around 2000 kms" 1 are consistent with our data. 
We also argue that the blending of the ubiquitous (in our 
data) feature at 2.26 /xm with the 2-0 R branch band head 
has to be taken into account when using this profile for 
velocity determinations. 

We note that comparing the spectra of the four other 
supernovae with CO emission with those of SN 1987A it 
is consistently the case that the spectral signature of the 
emission resembles that of SN 1987A at a later age than 
that at which the other SNe are observed. 

The other aspects of our data indicate that the spectra 
of type II SNe in the near IR are very similar and that 
recombination continues to play a major role even 5 to 6 
months after explosion. 

A better understanding of the physical conditions in 
the ejecta of supernovae can be achieved with higher qual- 
ity infrared data as would be obtained with ISAAC at the 
VLT. Very high quality data may even constrain the iso- 
topic ratios of carbon and oxygen as the wavelengths of 
the band heads to shift depending on the exact isotopes. 
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